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Abstract 


On  our  31st  trip  to  the  laser  facility  at  WSMR  we  carried  out 
experiments  on  laser  ablation  of  black  and  white  Delrin  [also  called 
polyoxymethylene,  polyformaldehyde,  (HCHO)x].  Mass  ablation  and 
thrust  generation  (Impulse)  were  accurately  measured  as  a  function 
of  input  laser  energy  in  one  shot  experiments.  The  efficiency  of 
conversion  of  laser  energy  to  jet  kinetic  energy  depended  on  the 
geometry  of  the  energy  absorption/conversion  zone.  The  most  ideal 
geometry,  an  axis  symmetric  mini  thruster,  produced  ~  60  % 
conversion  efficiency.  The  extensively  studied  10-cm  diameter 
Lightcraft  (with  inverted  paraboloid,  plug  nozzle  geometry)  produced 
~  50%  conversion  efficiency.  The  upper  limit  to  energy  conversion 
was  computed  with  CEA  code  to  be  73%  for  the  well  defined  mini 
thruster  geometry.  Thus,  total  losses  amount  to  ~  13%  and  ~  23%. 
This  is  a  significant  finding  and  helps  to  validate  the  concept  of 
“momentum  calorimetry”,  in  which  experiments  like  those 
accomplished  here  may  be  conducted  to  obtain  reliable  heats  of 
formation.  The  performance  of  candidate  chemically  enhanced  laser 
ablation  or  other  solid  propellants  may  be  measured  on  a  small  scale. 
In  these  most  recent  experiments,  a  near-exact  match  of  coupling 
coefficients  (1%)  was  achieved  in  a  14-fold  scale-down  of  the  10-cm 
Lightcraft  to  the  mini  thruster. 
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Outline 


*  Collaboration  Network/Why  Laser  Propulsion? 

*  Flavor  of  PLVTS 

*  The  Laser  and  EL  measurement  (Joules) 

*  The  Pendulum  and  I  measurement  (Newton  seconds) 

*  The  Mettler  Balance  and  m  measurement  (milligrams) 

*  Compare  the  EL,  I,  m  measurements  on  2  Test  Articles 
—  Light  Craft,  model  200-3/4 

—  Mini  Thruster  Standard  for  momentum  calorimetry/prop  devel 

*  Elementary  considerations  (energy/momentum) 

*  Comparison  of  experiments  to  1-D  equilib  code  (CEA) 

*  Conclusions/Work  in  progress/Flight  Tests  Movie 
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Phase  II  Program  Collaborations 
X-50LR:  Experimental  50-cm  Laser  Ramjet 


Team  December  2004 

Pulsed  Laser  Vulnerability  Test  System 


Frank  Mead,  Bill  Larson,  AFRL/PRSP 
Jim  Shryne,  RSI 

John  Harchanko,  POLARIS  TECHNOLOGIES 
Steve  Squires,  Chris  Beairsto,  Mike  Thurston 
JaySpray,  WSMR/HELSTF/PLVTS 


JOHN  P  STAPP 
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Overall  Energy  Conversion  in  Laser 

Propulsion  Mission 


Ef  =  l/2mv2  =  r|  a  (3  y  5  Ewall 


r|=  propulsion  efficiency  (jet  kinetic  energy  to  vehicle  kinetic  energy) 
a  =  expansion  efficiency  (internal  propellant  energy  to  jet  kinetic  energy) 

P  =  absorption  efficiency  (laser  energy  at  vehicle  to  internal  propellant  energy) 
y  =  transmission  efficiency  (laser  energy  at  ground  to  laser  energy  at  vehicle) 
8=  laser  efficiency  (electric  energy  to  laser  energy  at  ground) 

*****  Issue:  separability  of  rj  a  P  y  8  and  Ewall  ***** 

66  $500  worth  of  electricity  to  put  1  kg  into  LEO.” 

At  $0.1 0/KWH,  $500  buys  18,000MJ  =  Ewall  (1  KWH  =  3.6  MJ); 

1  kg  at  10  km/s  has  Ef  =  50  MJ,  so  rjaPyS  =  0.0028  =  50/18000 
But  if  28%  overall  efficiency,  then  $5/kg 


Phipps,  Reilly,  Campbell,  Laser  &  Particle  Beams  18  (2001)  661-695 
Pirri,  Monsler,  Nebolsine,  AIAA  Journal  12  (1914)  1254-1261 
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Pulsed  Laser  Vulnerability  Test  System 


Ifj  * 


Laser  Specifications 


Pulsed  C02  Laser 
10  KW 

~  5  to  30  ps  pulse  width 
Up  to  30  Hz 

Up  to  1000  J/pulse  (El  +  10%) 
Near  Field  Burn  Pattern 
~  1 0  feet 


Pulse  Shape 
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FFT  and  Far  Field  Burn  Patterns 


500  feet  1 000  feet  1 500  feet 
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10  cm 


Pendulum  Test  Stand 
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Ratio  of  Pendulum  Momentum  to  Hammer  Impulse 


Pendulum  Deflection  (degrees) 
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Hammer  Volts 


Traceable  Impulse  Calibration 


1275  1300  1325  1350 

time  (ms) 
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RVDT  Volts 


Mettler  Balance  or  Digital  Balance 


Measure  m  to  +  0.3  mg 
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Model  200  Lightcraft  Series: 

An  AF-Patented  Laser  Vehicle  Concept 
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and  Mini-Nozzle  Standard 


200-3/4  Lightcraft 


\/  10g  Delrin 

/  Parabolic 
v  Optic 

s(ideal  plug  nozzle)  =  14 
m=40g 

Delrin  surface  area  ~  25  cm2 
350  J/25  cm2/18  |us  =  0.8  MW/cm2 
Cm=450  N/MW,  EL/m=5.1  MJ/kg 
Ve=  2270  m/s,  efficiency  =  0.51 
T/W=CmP/mg  =  11  at  P=10  KW 


Mini-nozzle 
26°  divergence  angle 

5  cm 

- ► 

A 


s  =  16 


8=8 

m=7.8  g 

Delrin  surface  area  -  0.71  cm2 
25  J/0.71  cm2/1 8  ^is  =  2.0  MW/cm2 
Cm =442  N/MW,  EL/m=6.3  MJ/kg 
Ve=2795,  efficiency=0.62 
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Laser  Light  Craft  Flights 


+ 
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Light  Craft  Flights  with  Air  or  Delrin 
Model  200-3/4,  m=40g  lOkW  at  25  Hz  400  J/Pulse 
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Air  Plasma 


Distribution  A  -  Approved  for  public  release,  Distribution  Unlimited 


18 


Mini  Thruster  25  J,  18  |as,  0.71  cm2 


Black  Ser  9 
24.9  J,  18  (is  pulse 
3.8  mg/shot 


10.2  mNs 


10.9  mNs 


12.2  mNs 


11.1  mNs 


11.0  mNs 


White  Ser  10 
24.9  J,  18  (.is  pulse 
4.0  mg/shot 


AIR  Ser 11 
24.9  J,  18  (is  pulse 
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I,  m,  El  for  Mini  Thruster 


Miniature  Nozzle  w  ith  black  or  w  hite  Delrin 


laser  energy  (J) 


0  5  10  15  20  25  30 

laser  energy  (J) 


I/El  =  444  Ns/M  J  m/EL  =  0.1 60  mg/J 

Ve  =  (l/EL)/(m/EL)  =  2775  m/s 
Efficiency  =  y2(l/EL)2/(m/EL)  =  0.616  =  ap<P 
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10  cm  Light  Craft  322  J,  18  jas 


139.2  mNs 


144.8  mNs 


153.9  mNs 


146.7  mNs 


150.2  mNs 


154.1  mNs 


151.8  mNs 


133.0  mNs 


136.1  mNs 


137.8  mNs 


Black  Ser  7 
322  J 


White  Ser  6 
322  J 


59.8  mg/shot 


68.7  mg/shot 
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I,  m,  E  _  for  Light  Craft  200-3/4 


1 0  cm  w  hite  and  black  Delrin 


EL  (J) 


1 0  cm  w  hite  and  black  Delrin 


I/El  =  447  Ns/MJ  m/EL  =  0.201  mg/J 

Ve  =  (l/EL)/(m/EL)  =  2224  m/s 
Efficiency  =  y2(l/EL)2/(m/EL)  =  0.497  =  ocp<t> 
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CONVERSION  OF  LASER  ENERGY  TO 
JET  KINETIC  ENERGY 


Q*  = 


PEL 


in, 


Specific  internal  energy 


P 


Ejet  =  i  mp  <  ve2  >  =  ampQ*  =  a(5EL 
1=  mp<ve> 

-jF=ap<v*>2  =apo 

2mpEL  p  1  p 
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t  mf 

{Fdt  Jd(mve) 

0  
mf  mf 
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AIP  Conference  Proceedings, 
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Performance  map  of  known  laser 
materials  and  theoretical  air 


exit  velocity  vE  =  I/m  (m/s) 


Larson,  Mead,  Knecht,  AIAA  2004-0649 
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Experimental  Results  December  2004 
200-3/4  Light  Craft  and  Mini  Thruster 


0  1000  2000  3000  4000 

Ve(m/s) 
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apO  =  1/2CmV, 


I 


Obtain  P  -  20,000  bar  T  ~  3700  K  via  CEA 
Specify  expansion  ratio 
8  =  4,  8,16,  32,  64 

Obtain  Isp,  thermo  props  in  exit  plane  via  CEA 


Distribution  A- Approved  for  public  release,  Distribution  Unlimited 


26 


6  J/mg  Energy  Addition  to  Delrin 


HCHO  (g)  (u,p) - ►  CO(g)  +  H2(g)  (u,p)  +  1  mg  02 


420  s 
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8.8  J/mg 


lifetime  (s) 


Lifetime  of  Formaldehyde,  t(T,P) 


HCHO  +  M  =  H  +  HCO  +  M,  k=5x1015exp-308kJ/mol/RT  cm3/mol-s 
F.  Gernot,  D.  F.  Davidson,  R.  K.  Hanson,  Int.  J.  Chem.  Kinet.  36  (2004)  157 


1000 


Mechanism: 

(i)  HCHO  +  M  =  H  +  HCO  +  M 
(p)  H  +  HCHO  =  H2  +  HCO 
(t)  H  +  HCO  =  H2  +  CO 
(t)  H  +  H  +  M  =  H2  +  M 
(t)  HCO  +  HCO  =  CO  +  HCHO 


3000 


5000 


Temperature  (K) 
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Mole  Fractions  at  Equilibrium 

Formaldehyde  expansion  from  P=22694  bar,  T=3732K 
[rho=1227  kg/m3, u=1. 975  MJ/kg] 


Species 

chamber 

CO 

II 

CO 

s=64 

species 

chamber 

00 

II 

CO 

8=64 

CO 

0.47502 

0.48415 

0.35692 

CH30H 

0.00015 

0 

0 

H2 

0.39082 

0.39891 

0.36466 

CH3CHO,ethanal 

0.00014 

0 

0 

H20 

0.06058 

0.04282 

0.08215 

C3H4,allene 

0.00013 

0 

0 

CH4 

0.03818 

0.05811 

0.05318 

C3H6, propylene 

0.00013 

0 

0 

C02 

0.00856 

0.01574 

0.05707 

CH2 

0.00012 

0 

0 

C2H2 

0.00742 

0.00002 

0 

C2H2,vinylidene 

0.00009 

0 

0 

CH3 

0.00472 

0.00001 

0 

CH20H 

0.00007 

0 

0 

H 

0.00402 

0 

0 

C3H5,allyl 

0.00006 

0 

0 

C2H4 

0.00267 

0.00014 

0 

C4H2 

0.00006 

0 

0 

HCO 

0.00180 

0 

0 

COOH 

0.00005 

0 

0 

HCHO 

0.00180 

0.00003 

0 

CHCO.ketyl 

0.00004 

0 

0 

CH2CO 

0.00096 

0 

0 

CH30 

0.00003 

0 

0 

C3H3,2-pryl 

0.00039 

0 

0 

C2H 

0.00003 

0 

0 

C2H3, vinyl 

0.00035 

0 

0 

C302 

0.00003 

0 

0 

OH 

0.00032 

0 

0 

C4H6, butadiene 

0.00003 

0 

0 

C2H6 

0.00027 

0.00005 

0.00001 

C20 

0.00002 

0 

0 

HCOOH 

0.00026 

0 

0 

C2H50H 

0.00001 

0 

0 

C3H4 

0.00025 

0 

0 

C3H4,cyclo- 

0.00001 

0 

0 

CH3CO 

0.00019 

0 

0 

C3H8 

0.00001 

0 

0 

C2H5 

0.00019 

0 

0 

C4H6,1butyne 

0.00001 

0 

0 

C(gr) 

0 

0 

0.08601 
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Mole  Fractions  at  Equilibrium 

Formaldehyde  expansion  from  P=230  bar,  T=3433  K, 
[rho=12.02  kg/m3,u=1.907MJ/kg] 


mole  fractions 

Chambr 

throat 

8=4 

e=8 

e=16 

e=32 

e=64 

CO 

0.49263 

0.49587 

0.49955 

0.49553 

0.47166 

0.43233 

0.38881 

H2 

0.47969 

0.48865 

0.49681 

0.4913 

0.48006 

0.46937 

0.46168 

H 

0.02313 

0.01223 

0.00003 

0 

0 

0 

0 

H20 

0.00239 

0.00173 

0.0014 

0.00435 

0.01321 

0.0231 

0.03105 

C2H2, acetylene 

0.00105 

0.00074 

0.00002 

0 

0 

0 

0 

CH4 

0.00032 

0.00033 

0.00176 

0.00435 

0.00673 

0.00753 

0.00728 

C02 

0.0003 

0.00023 

0.00041 

0.00224 

0.01093 

0.02605 

0.04371 

CH3 

0.00021 

0.00012 

0 

0 

0 

0 

0 

HCO 

0.00014 

0.00005 

0 

0 

0 

0 

0 

*OH 

0.00005 

0.00002 

0 

0 

0 

0 

0 

CH2 

0.00002 

0.00001 

0 

0 

0 

0 

0 

HCHO,formaldehy 

0.00002 

0.00001 

0 

0 

0 

0 

0 

C2H 

0.00001 

0 

0 

0 

0 

0 

0 

C2H2,vinylidene 

0.00001 

0 

0 

0 

0 

0 

0 

C2H4 

0.00001 

0.00001 

0 

0 

0 

0 

0 

C(gr) 

0 

0 

0 

0.00223 

0.01741 

0.04162 

0.06748 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (e) 


Pc=20.2  bar,  Tc=3060  K,  pc  =1.18  kg/m3 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (s) 


expansion  ratio 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (e) 


e  expansion  ratio 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (e) 


s  expansion  ratio 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (e) 


s  expansion  ratio 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (e) 


e  expansion  ratio 
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Blowdown  from  specified  initial  state 
(u,p)  with  specified  expansion  ratio  (e) 


CO 


s  expansion  ratio 
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Experimental  data  (I,  EL,  m)  and  derived 
parameters  (Cm,  Ve,  efficiency,  EL/m) 


Geometry 

I  vs  El 

slope 

m  vs  El 

slope 

Cm 

Ve 

Efficiency 

EL/m 

mNs/J 

R2 

mg/J 

R2 

Ns/MJ 

m/s 

MJ/kg 

Mini  thruster 
white 

0.444 

0.97 

0.160 

0.99 

444 

2775 

0.616 

6.3 

Mini  thruster 
black 

0.439 

0.98 

0.156 

0.98 

439 

2814 

0.618 

6.4 

Mini  thruster 
AIR 

0.253 

0.97 

- 

- 

253 

- 

- 

- 

10-cm  Model 
white 

0.447 

0.85 

0.201 

0.96 

447 

2224 

0.497 

5.0 

10-cm  Model 
black 

0.453 

0.92 

0.194 

0.93 

453 

2335 

0.529 

5.2 
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Conclusions/Work  in  Progress 


•  Cm=450  N/MW  for  Light  Craft/Delrin  (350  J,  18  (as) 

•  Cm=442  N/MW  for  Mini  Thruster/Delrin  (25  J,  18  (as) 

•  51  %  efficiency  for  EL  to  jet  KE  for  Light  Craft 

•  62  %  efficiency  for  EL  to  jet  KE  for  Mini  Thruster 

•  Future  Experiments 

—  Vary  pulse  width,  5  and  30  jus,  expansion  ratio,  8  =  4,  16,  ... 

—  Increase  EL  up  to  -  100  J/pulse  in  mini  thruster 
—  Measure  time  resolved  thrust  with  piezoelectric 
—  Develop  chemically  enhanced  ablative  propellants 

•  Future  Calculations  with  Chemical  Equilibrium  Applications  Code 

—  Factor  pressure  thrust  into  analysis 
—  Analyze  Chemically  Energetic  Propellants 

Distribution  A- Approved  for  public  release,  Distribution  Unlimited 


39 


O  for  Bimodal  velocity  distribution 


Chunks  of  propellant 
Hot  gases 


fheavy  mass  fraction,  vslow  velocity 
f,ight  mass  fraction,  vfast  velocity 


<V>2  =  (fheavyVsIow  +  fligh.Vfast)2 

®  =  <V>2/<V2>  =  (fheavy  +  flightr)2/(f. 


<V2>  =  f, 


heavy  v  slow 


2  +  f, 


lightvfast 


+  fiiah.r2)  where  r  =  vfast/vslow  >  1 


heavy  light 


r  =  high  velocity /low  velocity 


r  =  600 


mass  fraction  of  heavy  m  ass 


r  =  4000 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 
mass  fraction  of  heavy  mass 
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Ratio  2D  or  3D  to  ID 


Divergence  Loss 


i 


0.98 


0.96 


0.94 


0.92 


□  NASA  1960 

— a — 3D:  0.5(1  +cos 

— h — 2D:  (sin  a)/a 

\c 

3V 

a) 

LJ  1  1 

1  LJ  1 

J 1  LJ 

LJ  1  J 

1  LJ  1 

J 1  LJ 

0  5  10  15  20  25  30 

a  =  V2  (divergence  angle) 
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